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Abstract 

The shape of drop on a flat horizontal plane is obtained by including 
the first order of correction by the weight. The sphere solution of the 
weightless drop is used to introduce a new polar coordinate by which the 
perturbative expression for a region of a drop can be extended analytically 
to the entire surface of a drop having both the concave and the convex 
parts. Comparison with experimental data are presented. 
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1 Introduction 



To the interface of two mediums a and b, it is assigned an energy per area 
of interface, the so-called interfacial energy coefficient 7 a b. For example, the 
liquid-vapor parameter 7i v = 7 describes the energy content coming from the 
fact that liquid molecules near the surface have less neighbors than those in the 
bulk. The corresponding coefficient is called the surface tension. The shape of 
a drop of liquid on a solid surface, in the idealized case (absence of impurities 
and pinning effects), is determined by the quantities: 1) the surface tension 7, 
2) the adhesion coefficient er, 3) the shape of the solid surface, and due to the 
weight, 4) the drop's volume. The adhesion coefficient is defined by the surface 
tension, and the solid-liquid and the solid- vapor interfacial energies as [TJ [2] 

er = 7 S v + 7 - 7si- (1) 

At the solid-liquid- vapor point of contact, the contact angle $ in the equilibrium 
condition is given by the Young equation 

cosi?=--l. (2) 

7 

The above relation determines 1? for < a < 27. Three classes of possibilities 
for the contact angle are presented in Fig. 1. The cases with a « 27 and a <C 7 
correspond to the highest and the lowest spreading of the drop on the solid 
surface, respectively. Hence the name complete wetting befits to the case with 
a k, 2 7 . 
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Figure 1: Three classes of the drop's shape on a solid surface. 



At every point on the drop surface the Young-Laplace relation holds [TJ |2] 

in which Ap = p\ — p v is the difference pressure across the surface, and (Ri, R2) 
are two principal radii of curvature of the surface at the point. At each point 
of the drop's surface the total curvature R^ 1 + R^ 1 is determined in terms of 
the surface equation and its derivatives. Provided by the hydrostatic laws, Ap 
can be expressed in terms of the surface equation as well. So the Young-Laplace 
relation is the partial differential equation which, accompanied by appropriate 
boundary conditions, determines the shape of the drop's surface. 

Equivalently, the shape of drop can be obtained by the minimizing the energy 
of a static system. While the surface tension tends to decrease the surface area 



f 



of the drop, the adhesion coefficient tends to increase the surface area of the 
contact region, and the gravity tends to lower the center of mass of the drop. The 
competition between these effects determines the shape of the drop. For a drop 
with volume V and density g, one can give an estimation for each of the effects. 
The order of the drop's size is estimated by L — V 1 ^ 3 . In many practical cases 
the surface tension and the adhesion coefficient, though with opposite effects, 
may be considered at the same order, meaning 7 and a are comparable. So 
the contribution of the interfacial energies, which is proportional to the area, 
is estimated by 7 L 2 . The contribution of the weight to the energy is given by 
gVgL = ggL 4 , with g as the gravitational acceleration constant. Comparing 
these contributions, one can differentiate three regions: 

• L <C \/j/gg- the effect of weight is small; 

• L ~ \flTog~'- the weight and the interfacial energies have comparable 
effects; 

• I> \Jll gg~- the effect of weight is dominant. 

In other words, the comparison between the length i = y / t / ' gg and L, or equiv- 
alently the value of the Bond dimensionless parameter L 2 gg/j would determine 
the regime. For a weightless drop only the contribution from the surface tension 
exists. Minimizing the area for a fixed volume, the shape of the drop's surface 
turns out to be part of a sphere. 

The problem of a drop on a horizontal surface with the effect of surface ten- 
sion being balanced with gravity has been studied for more than a century. The 
early numerical solutions go back to 1883 [3], with updates by different authors 
[H [51 [5] . Different perturbative treatments of the problem have been developed 
over the years, among them are those by El El [XD] for small drops (small Bond 
number). As large drops (or vanishing surface tension drops) are theoretically 
an infinitely large and thin film of liquid subjected to the boundary conditions 
at the outer edge, the limit of large Bond number falls and has been studied in 
the context of singular perturbation problems [11] . Based on the similarity be- 
tween the truncated oblate spheroid and drop's shape, an approximated profile 
is suggested in ;12] for the shape of the drop. In |13] a new numerical treatment 
of the problem is given based on a variational method to minimize the total 
energy of the drop, by which the use of the tables by [3] is more direct than 
the earlier treatments. As another effort in this direction, in [T5] the singular 
perturbation technique is used to obtain the asymptotic expressions describing 
the shape of small sessile and pendant drops. The study of the profiles of resting 
drops in different situations are particularly important for practical purposes. In 
fact, one of the most common methods to measure the surface tension of liquids 
is based on the matching between calculated drop's profiles and the measured 
drop's shapes. Over the years, the optimization of matching methods between 
the calculated profiles and the experimental data on drop's profiles has been the 
subject of several research pieces P~5l fTBl fTTj . 
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It is the purpose of this note to develop a perturbation method to include 
the effect of the gravity on the shape of the drop's surface. Here in particular 
the shape of drop is obtained by including the first order of correction by the 
weight. The sphere solution of weightless drop is used to introduce a new polar 
coordinate by which the perturbative expression for a region of a drop can be 
extended analytically to the entire surface of a drop having both the concave 
and the convex parts. 



2 The shape of weightless drop 

It is known that the shape of a weightless drop is a part of a sphere. A proof 
for this is, however, given here, mainly to introduce the mathematical tools for 
the perturbative method. In particular, it is the sphere solution by which the 
formerly mentioned polar coordinate is defined. 

It can be shown that the total curvature of the surface z = f(x,y) is given 
by [IB] 

1 1 (1 + f 2 y )h X + (1 + fDfyy - HJyhy 

Ri R2 (i + 11 + flf' 2 ' 1 ' 

in which f x = d x f, f y = d y f, f xx = d^f, f yy = d y f and f xy = f yx = 
d x d y f. In the case with cylindrical symmetry, / depends on the combination 
p = \J x 1 + y 2 . In this case we find 



1 , lid/ |f 



Ri R 2 p dp Y ^l + f' 2 J ' 



(5) 



where /' = — . Defining 
dp 



Hp) = , f (6) 



f" 

by which J> (p) = (1 + / , 2)3/2 , we get: 

7r + 7r = i -f O^D- (7) 

Ri R 2 pdp K 

In absence of gravity (g = 0), using Q, we find for the Young-Laplace relation 

in which Apo is a constant, representing the difference pressure due to the surface 
tension of liquid in absence of gravity. In fact, due to the surface tension, the 
pressure inside the liquid drop is bigger than outside. So Apg, and thus kq, are 
positive. 
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Figure 2: The geometry of the mathematical setup for ■& > 90°. 



In the case in which both concave and convex parts are present the entire 
surface can not be represented by only one function depending on p, simply 
because there are two different values of z for some p (Fig. 2). So the equation 
([8| should be solved for the convex (/o+) and concave (/o-) parts separately 
(Fig. 2). We mention |/q ± | = T/o±- I n the zero gravity case, however, since 
the right-hand side of equation Q is constant, one solution is related to the 
other simply by changing the sign of its derivative, or its corresponding tpo by 
((6). With the integration of (|8) 

T pipo± = K Q p 2 + a± (9) 

where a± are the constants of integration. By the above reasoning, in the present 
case a + = a_. Also a + should be set to zero in order that -00+ does not blow 
up at p — 0, for which by definition IVKp)! < 1. One then has 



Jo 



12 

± _ .,2 2 



and so 



<P 2 (10) 



dp ^l^nfp- 
for which by the integration we find 



d/0± = T^J^ (ID 




(12) 



which represents a sphere with radius R = k^ 1 > whose center is located on 
z-axis at z — zq. As mentioned, in the case with i? < 90° only the positive sign 
has meaning, and both ± signs should be kept in the case with d > 90° . 

It is useful to check for the number of parameters involved. The parameters 
Apo, or equivalently Ko, and zo are unknown at the first place. Following a 
simple geometrical argument in the sphere (see Fig. 2), we have 

co S1 ? = -|, (13) 
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by which we have z > for i? > 90° and z < for d < 90°, corresponding to 
center above and below the solid surface, respectively. So given by the relation 
for the volume 

V = -R 3 (l- cosi?) 2 (2 + cost?) (14) 
3 

zq and R are fixed. The parameters po and p\ in Fig. 2, as the contact and 
equatorial radii, respectively, can be obtained once the equation of sphere being 
solved for z = and z = z , yielding 



Po = iZsim?, pi = R- 

The place of drop's apex in the spherical solution is 

h Q = R + z Q 



(15) 



(16) 



Before to proceed, let us to introduce an identity for later uses. The dif- 
ference pressure in presence of gravity gets contribution from the weight of the 
drop's layers as well. So we have for the ratio 



Ap(z) 

7 



2k 



7 



(h-M) 



(17) 



in which h is the height of the drop's apex, and k, similar to its counterpart k , 
represents the difference pressure due to the surface tension but here in presence 
of gravity. The contribution of the surface tension to pressure in presence of 
gravity is different and simply can be understood as the area of the drop is 
changed due to the effect of weight. Hence the surface tension contribution to 
the energy of the drop is different due to the gravity. So, the Young-Laplace 
relation reads 

T~(ptf±)=2K+^(fc-/±). (18) 
p dp 7 

Integrating the above for the upper and lower parts gives the followings, 

rPi 



Pi 

pi - po sin <Q 
in which we have used 



27 
Q9_ } 

27 



2 Q9 

Pi , 

7 Jo 

(pI-pI) 



Pf+(p)dp 
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1 



pf-(p)dp 



f- 



1 + f'l 



— tan$ 
v 7 ! + tan 2 d 



= sin?? 



(19) 
(20) 

(21) 



for ■d > 90°. Subtracting g9j and @ gives 

sint? 



01} 

27 



Po 



ggv 

2tt7Po 



(22) 
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in which we have used the relation for the volume of drop, 



V 
2^ 



Pi 



P/+(p)dp 



Pi 



Pf-(p)dp 



(23) 



/Hi 



It is easy to show that identity ( 22 ) is valid for the acute contact angle (# < 90°), 
as well. It is reminded that in obtaining ( 22 1 no approximation is used, and so 
it is an exact relation. 



3 The shape of lightweight drop 



Here we consider the first correction of gravity to the shape of a drop, supposedly 
applicable to the drops with tiny weight, or equivalently small volume. It is 



clear from (181 that the effect of weight, as mentioned earlier, appears in the 
combination gg/j. By the help of volume V, it is useful to introduce the Bond 
dimensionless parameter A = V 2 / 3 gg/~/. This supposedly small parameter helps 
to develop a perturbative expansion for the contribution of gravity on the shape 
of drop. At the first order of correction one has 

z = f(p)=f (p) + Xf 1 (p) (24) 

where fo(fl) is the sphere solution found in the previous section. By inserting 
above in ([6| one finds for ip+(p) 

f+ fo+ , x fi+ 



-0(A 2 ). (25) 



The combination k + in the right-hand side of ( |18[) can easily be rearranged 
based on the effect of gravity using the identity ( |22p . As it is expected that 
the contact radius po is changed under the effect of gravity, at the first order 
of perturbation in g, by replacing po = Rsm.'d + Spo, the identity comes to the 
form 



27 



1 

R 



Spo 
R 2 sintf 



ggV 

27T 7 R 2 sin 2 



(26) 



We later will find that 8p > 0, as expected. By insertion (25) and (26) in (18), 
and using the fact that /o satisfies the equation with A = (eq. (8f), for the 
case with -d < 90° or the upper half of case with d > 90° the Young-Laplace 
relation reads 



d 

dp 



1 



in which 



(l + /&) 3/ 
27 Spo 



1/2/3 



P{a + fo+{p)) 



V 



gg R 2 sin $ irR 2 sin 2 # 



Using f 0+ (p) = ^R 2 - p 2 



1 



( 1 + /o+) 3/2 1^ 2/3 V2 



zq, integrating (27) from to p gives 
(\{z* + a)p 2 -\{R 2 -p 2 r 2 + \R3 



(27) 



(28) 



(29) 
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Again using the expression for /q+ (p) , we find 



fi+(p) 



R 3 

y2/3 



(z + a )p 
2(R 2 - p 



2)3/2 



1 

3p 



R 3 

3p{R 2 ~ p 2 ) 3 / 2 



(30) 



It is useful to note that for the above expression, the limit p — > exists and is 
zero, as expected. However, we mention that the above expression diverges in 
the limit p — > i?, indicating that the present form of the perturbative solution 
fails for drops with contact angle close to or greater than 90° . We will come back 
to this issue later. The above expression can be used to find the corrected value 
of the contact radius po, or equivalently 6po- Subjected to the condition that the 
drop intercepts the solid surface with contact angle and also the constraint 



on the volume of drop, the expression (30 1 should satisfy the fallowings 



tantf = f'(p ) ~ /' (p ) + \f'(Rsm§), 



(31) 



V = TT 



Po rPo /-ilsintf 

p 2 \f' + (p)\dp^-7r / f?fL(f)dp-wX / P 2 f 1+ (p)dp (32) 



(i 



i) 



in which po = R sini? + Spo. We mention, due to presence of A, it is sufficient to 
insert the unperturbed values of previous section in the expression for f[ + . It is 



easy to check that, thanks to the identity (22 ), the first in above is automatically 



satisfied. By the second condition the change in the contact radius is found to 
be 

ggR 3 (1-cos^) 2 
SP0 = 6 7 sintf(2 + costf) >0 (33) 

by which the constant a is obtained 

,(l-cosi?)(3 + cosi?) 



-R- 



3(2 + cost?) 



(34) 



Once again the relation (30) can be integrated, leading to 



h+{p) 



R 3 

V 2/3 



3{z 



2R 



1 (R 
3 ln! 



Vr^ 



2R 



+ b 



(35) 



where b is a constant, that should be determined by the condition 



= f+(po) ~ / 0+ (po) + A/ 1+ (i?sin^). 



(36) 



Using above we find 



b = 



R 3 



V 2/3 



cos $ 



2(2 + cos ■#) 



2 
-Incos- 



(37) 



In the case i9 > 90 

of the drop. It is seen that both (j30h and (|35|) diverge as p 



( 35 1 represents the correction only for the upper half 

R. So the 



perturbative solution in the present form fails at the vicinity of the circle p = R. 
For d > 90° one has to find the solution for the convex part (lower half part) 
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z 




Figure 3: The geometry with the polar coordinate. 



as well, represented by the f\~- This also can be done in similar lines of the 
concave part, however the boundary conditions are different, ft is easy to see 
that in this case, just like we had for the upper part, fi-(p) diverges for p — > R. 

One way out of the divergent behavior of fi±{p) near p = R is to change 
the role of function z and the variable p, working with p = h(z), for p > 0. One 
can then find valid perturbations near p = R (the maximum of p), however, 
this time the solution diverges at the top of drop (p — > 0). This shows that 
the appearance of the divergent behavior in the perturbative expressions for the 
shape of drop is not an intrinsic one, but a coordinate artifact. By these all, 
there are three functions /i±(p) and hi(z), which should be joined smoothly to 
give the correction of the gravity to the shape of the drop in all regions. 

A better way is to use the polar coordinate suggested by the sphere solution 
of the weightless drop. Choosing the bottom of the circle corresponding to the 
weightless drop as the origin, and measuring the angle from the z-axis, one has 
(Fig. 3) 



r(6) cos9 - d = z = f( P ) = f (p) + Xfi(p) 



(38) 

where d = R — z = R(l + cos??), by |T3"1 ). A perturbative expansion for r(9) is 
then 

r{6) = r (6) + \n(6) (39) 



where 
Putting 



r-o (0) = 2Rcos6. 
p = r(6) sin = R sin(20) + \r x (0) sin 



in ( 38 1 , one obtains 



A : r o (0)cos0- J R. + z o -/o(p) 
A 1 : ri(0)cos0 = 



p=Rsin(20) 



ri (0) sin 0^ +/l ( p ) 



(40) 
(41) 

(42) 
(43) 



p=_Rsin(29) 
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As the concave part always exists, we can use /o+ and /i+ in above to find the 
unknown part J~i(0), 



n(0) 



R 3 
V 2 / 3 cos( 



3(z 



ha) 
6^ 



2R 



cos(20) In cos0 



cos(20) , . 
cost/ 



in which constant a and b are given by (34 1 and (371, respectively. In fact the 



above expression is nothing but the analytically extended result for the concave 



part ( 35 ) to the entire surface of the drop. We mention ( 44 1 has smooth behavior 



for the whole interval < 9 < w, which covers the convex part too. It is useful 
to define the angle 9o, as the polar angle at which the contact of the drop and 
the surface takes place. This angle is easily obtained by the following condition: 



d tan 9q = po ~ i?sin$ + Spo, 



(45) 



in which Spo is given by (33). For spherical solution 8q = $/2. 



It is a matter of interest to obtain the equatorial radius (maximum bulge) 
Px for the case with d > 90°. It is obvious that by spherical solution pi = R, 
happening at the angle 9\ = 45°. In general, the equator is defined by the 
condition: 



= p'(0i) = 2 J Rcos(26»i) + A 



AO 



(sin0ri(0)) 



0=?r/4 



Using (41 ), at the first order of A, we find 

Pi 



IT /IT 

R + Asm — n - 
4 \4 



The equatorial plane intercepts z-axis at 
Z\ = r{6\) cos Ox — d 



zo 



i?cos(20i) + A cos - ri^-j) 



(46) 

(47) 

(48) 
(49) 



which is easy to find by using (46). For later use, the distance between the apex 



of drop and the equatorial plane, h, is given explicitly 



h = h — z\ = R 




(50) 



4 Comparison with data 

In order to asses the accuracy of the perturbative solution presented in this 
work, here the outputs of the solution are compared with some available data. 
However, it would be helpful to begin with a demonstration of the results. 
In Fig. 4 the sphere and the perturbative solution are plotted for two drops 
which are equal except in their contact angles. As expected, the apexes of the 
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a) b) 

Figure 4: Two demonstrations of the sphere solution (dashed line) and the 
perturbative solution (solid line) for drop a) with acute contact angle (# = 
80°), and b) with obtuse one (0 = 125°). For both drops: V = 0.025 cm 3 , 
q = 1 g/cm 3 , 7 = 70 dyn/cm, g = 980 cm/s 2 . 



sample 


vol. 


cont. 

angle 


experiment 
PO 21 


sphere 

Po 
(err.) 


sol. (l2} 

Zl 
(err.) 


perturb. 

PO 
(err.) 


sol. pi} 

(err.) 




units 


















1 


72.0° 


0.1748 


0.1148 


0.1713 


0.1245 


0.1742 


0.1196 


water on 










(2.0%) 


(8.4%) 


(0.3%) 


(4.1%) 


carbon 


2 


71.3° 


0.2240 


0.1411 


0.2171 


0.1557 


0.2229 


0.1460 


steel [12] 










(3.1%) 


(10%) 


(0.5%) 


(3.4%) 




3 


71.2° 


0.2360 


0.1565 


0.2487 


0.1780 


0.2575 


0.1634 












(5.4%) 


(14%) 


(9.1%) 


(4.4%) 


water on 
PMMA [T71 


cm'* 
0.1234 


73.44° 


0.4897 




0.4462 
(8.9%) 


0.5007 

(2.3%) 





Table 1: The experimental data and theoretical values by spherical solution and the 
perturbative one for drops of water on two different surfaces. For carbon steel data: 
1 unit of vol.=6.75 x 10~ 3 cm 3 , 7 = 72 dyn/cm. For PMMA data: 7 = 70.6 dyn/cm. 
For both surfaces: g — 1 g/cm 3 , g = 980.7 cm/s 2 . All lengths are in cm. 



perturbative solution are lower than the sphere one's, while the contact radii 
are increased in comparison with the sphere solution. Also for the case with 
■d > 90° , the equatorial radius is larger than the radius of the sphere solution. 

The outputs of the perturbative solution are compared with data for drops 
with both acute and obtuse contact angles. In Table 1 the collection for the 
case with $ < 90° for water drops are presented. The data include three drops 
on Carbon Steel surface by [12 , and one drop on poly(methyl methacrylate) 
(PMMA) surface by [T7j . The comparison between the experimental values and 
the ones by the perturbative solution, perhaps except for the contact radius of 
third sample on Carbon Steel, shows a satisfactory agreement. The failure of a 
good agreement with the third drop might be better justified once the peculiar 
behavior of data with this drop is mentioned. For example, the contact radius 
of this drop is even less than the one with sphere solution, and the height of 
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10 3 x 


cont. 




experiment 




sphere sol. 


m ~ 


perturb, sol. 


|39| 


vol. 


angle 


in 


Po 


h 


Pi 


Po 


h 


Pi 


Po 


h 


cm 3 










(err.) 


(err.) 


(err.) 


(err.) 


(err.) 


(err.) 


0.370 


131.1° 


0.0445 


0.0337 


0.0442 


0.0458 


0.0345 


0.0458 


0.0462 


0.0357 


0.0456 












(2.8%) 


(2.3%) 


(3.5%) 


(3.9%) 


(5.9%) 


(3.2%) 


2.510 


129.5° 


0.0907 


0.0722 


0.0813 


0.0869 


0.0671 


0.0869 


0.0902 


0.0748 


0.0860 












(4.2%) 


(7.1%) 


(6.9%) 


(0.6%) 


(3.6%) 


(5.7%) 


4.830 


132.6° 


0.1163 


0.0884 


0.1035 


0.1074 


0.0791 


0.1074 


0.1137 


0.0957 


0.1057 












(7.6%) 


(11%) 


(3.8%) 


(2.2%) 


(8.3%) 


(2.1%) 


10.370 


132.4° 


0.1536 


0.1191 


0.1299 


0.1386 


0.1024 


0.1386 


0.1521 


0.1379 


0.1349 












(9.7%) 


(14%) 


(6.7%) 


(1.0%) 


(15.8%) 


(3.9%) 



Table 2: The experimental data and theoretical values by the spherical solution and 
the perturbative one for drops of mercury on glass slide by [Sj. q = 13.55 g/cm 3 , 
7 = 476 dyn/cm, g = 980.7 cm/s 2 . All lengths are in cm. 



apex is larger than the one of sphere solution. Both of these observations are 
in opposite with expectations from the effect of gravity on a real drop. 

In Table 2 the collection for the case with ■d > 90° for mercury drops are 
presented. The data include four drops of mercury on glass slide by [Sj. Here 
the data are available for the equatorial and contact radii, and the equator-apex 
difference h of the previous section. Again, perhaps except for the first drop, 
the agreement between experimental data and theoretical results is saisfactory. 
Just like the third water drop, also for the first drop here the data is in opposite 
expectations with the effect of gravity (less radii and larger height than the 
sphere solution) . It might be notable that in this collection the Bond parameter, 
as the perturbative expansion parameter, exceeds one (A = 1.33) for the last 
drop as the largest one. 
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